u 

< 
O 



> 

in 



THE DETERMINATION OF THE 'DIFFUSION COEFFICIENTS' AND 
THE STELLAR WIND VELOCITIES FOR X-RAY BINARIES. 
V.M.Lipunov 1 > 2 and S.B.Popov 1 
1 Department of Physics, Moscow University, 119899, Moscow, Russia, 
- 2 Sternberg Astronomical Institute, Universitetsky pr., 13, Moscow, Russia 

Q\ • e-mail: lipunov@sai.msu.su and popov@sai.msu.su 

G\ ' 

Abstract 

The distribution of neutron stars (NS's) is determined by sta- 
tionary solution of the Fokker-Planck equation. In this work using 
the observed period changes for four systems: Vela X— 1, GX 301-2, 
Her X-l and Cen X-3 we determined D, the 'diffusion coefficient ',- 
parameter from the Fokker-Planck equation. Using strong depen- 
dence of D on the velocity for Vela X-l and GX 301-2, systems ac- 
creting from a stellar wind, we determined the stellar wind veloc- 
. ity For different assumptions for a turbulent velocity we obtained 

V = (660 — 1440) km ■ s _1 . It is in good agreement with the stellar 
q ■ wind velocity determined by other methods. 

I/"") . We also determined the specific characteristic time scales for the 

'diffusion processes' in X-ray pulsars. It is of order of 200 sec for 
i-£h ■ wind-fed pulsars and 1000-10000 sec for the disk accreting systems. 

q ■ Keywords: accretiomneutron stars-stars:stellar wind-stars. 

CO 

1 Introduction 

> ■ 

• i-H , 

The most precisely determined characteristic for accreting neutron stars 
(NS's) is their period. Thus using observations of the period we can de- 
termine different properties of the observed object. 

Period changes show fluctuations. These fluctuations were discussed in 
de Kool Sz Anzer (1993). The authors determined noise level in these systems 
and characteristic time scales. Using these results we can estimate diffusion 
coefficients in the Fokker-Planck equation (see below) and the stellar wind 
velocities (only for the wind-accreting systems ). 

During accretion the angular momentum of plasma is transfered to the 
NS. But the process of the momentum transfer is not stationary. The trans- 
fered angular momentum fluctuates and therefore the period changes of the 
NS will also show fluctuations. 
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Table 1: 





p, sec 


Por&j sec 


L, erg/sec 


/i/10 30 Gs-cm 3 


^suobservi Y^S 




Vela X-l 


283 


7.7 • 10 5 


1.5 • 10 36 


3 


3000 


3000 


GX 301-2 


696 


3.6 • 10 6 


10 37 


120 


> 100 


100 


Her X-l 


1.24 


1.5 • 10 5 


10 37 


0.6 


3 • 10 5 


8000 


Cen X-3 


4.84 


1.8 • 10 5 


5 • 10 37 


5.7 


3400 


600 



Processes with fluctuations are well known (see for example Haken (1978)). 
Some applications of stochastic processes in astrophysics especially in accret- 
ing systems were discussed in Lipunov (1987),Hoshino & Takeshima (1993) 
and Lipunov (1992). In Hoshino & Takeshima (1993) the authors, using sim- 
ple models of MHD turbulence, try to explain aperiodic changes in X-ray 
luminosity of X-ray pulsars. Luminosity fluctuations are explained as the 
result of density fluctuations due to turbulence in the plasma flow. The 
authors used 2D model for accretion disk and 3d model for the wind accret- 
ing systems. Detailed exploration of this question is very difficult in both: 
theoretical and observational ways (there is no good theory of MHD turbu- 
lence and the resolution of modern equipment of satellites is not high enough 
for power spectra of X-ray pulsars (see Hoshina & Takeshima 1993). But 
detailed exploration of the density fluctuations will help to understand pe- 
riod fluctuations. It will be very interesting to compare X-ray luminosity 
fluctuations with fluctuations of the period of the NS. 

There are different methods of describing of these processes. In this work 
we use differential equations for the distribution function. For the frequency 
changes we can write the Langeven equation, which describes the process 
with fluctuations: 

^ = F{») + *(t) (1) 

Here, F(u>) -constant angular momentum. For F(u>) in the most general 
form we can write (Lipunov 1982): 

F(\omega )\cdot I=\left\{ 
\begin{ array} 

& \dot{M}\eta _k \0mega R_G"2-k_t\f rac {\mu ~2}{R_c~3}, 
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\qquad {\text wind\, accretion}-\\ 

& \dot{M}\sqrt{GMR_d}-k_t\frac {\mu "2}{R_c"3}, 

\qquad {\text accretion\, disk},\\ 

\end{ array} 

\right . 

where -R^-disk radius, k t and r^fc-dimensionless parameters {k t ~ 1, r/k ~ 1), 
/-moment of inertia of the NS, M- mass of the NS, fi-orbital frequency of 
the system, R G - radius of gravitational capture {Rq = J+^z , we include 

orb 

1 /3 

here the orbital velocity) and i? c -corotational radius, R c = (jjjr^ 

We assume that the 'force' is conservative and in this case we can write 

F(uj) in the form: F(u) = — V^V , where V is a scalar potential. $ is a 

fluctuating moment, i.e. < $(t) >= , (Lipunov 1987). 

The distribution of frequency, cu , is described by the distribution function 

f(co). This function satisfies the Fokker-Planck equation (Haken 1978): 

df _ dfF{u) | D d 2 f^ 

dt doj duj 2 ' 

where D is the 'diffusion coefficient', which is determined by the correlation 
of the stochastic force $ : 

<&(t)$(t')>=2D5(t-t') (3) 
Stationary solution of the eq.(2) is the following: 

f(u)=Nexp(-V(u)/D) (4) 
where N is determined from the normalisation condition: 

f(u)du=l (5) 

Using expression for V(uj) from Lipunov (1992) we can write D in the 
form: 

D = J±vL.^ (6) 

3GMI 7 V ; 
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Vela X-l 
GX 301-2 



A 

-9.1 
-8.5 



10 36.8 
10 37 



Table 2: 



where /c t -constant parameter {k t ~ 1) , /-moment of inertia of the NS and 7 
is evaluated as 7 ~ % , here t sn -time of spin-up and At -the characteristic 
time for period changes (see for details Lipunov (1987) or Lipunov (1992)). 

Using eq.(6) in 2.1 we shall determine the value of the 'diffusion coeffi- 
cient', D. With these D in 2.2 we shall make the estimates of the stellar wind 
velocities for Vela X-l and GX 301-2. 



2 Results. 

2.1 Determination of the 'diffusion coefficient'. 

At first we shall estimate D, using eq. (6). In this equation all variables, 
except At , are known (in principle). Their values taken from Lipunov (1992) 
are shown in table 1. 

Nagase (1992) gives results of GINGA observations of the cyclotron lines 
in X-ray pulsars. For Vela X-l and Her X-l values of magnetic field, B, 
coincides with the values \i that we used when radii of NS's are 10 km and 
6 km correspondently. The value for B obtained from observations for Vela 
X-l, B = 2.3 • 10 12 Gs, coincides quite well with the assumption that there 
is no stable accretion disk in this pulsar. 

Characteristic time At for the wind-accreting systems can be determined 
from the equation: 

At « 1.7 • 10 4 a- 2 10 2(A+8 - 5) L^fiJ sec, (7) 

where a - a fraction of the specific angular momentum of the Kepler orbit 
at the magnitospheric radius, A-noise level (see table 2) (de Kool & Anzer 
1993). 

In eq.(6) for t su we must use minimum values. These times were calculated 
using equations for pure spin- up from Lipunov (1992), t sumin are shown in 
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table 1. 

We took k t — 1/3, M — 1.5 M , / = 10 45 g ■ cm 2 . From eq. (7) we can get 
At for Vela X-l and GX 301-2. For Her X-l and Cen X-3, At is determined 
from the graph in de Kool & Anzer (1993) (see table 3). So we can write 
equation for D in the form: 

D = 5.55 • lO-^V^ 1 (yj^) " s ~ 3 ■ (8) 

Values of D for four systems are shown in table 3. 
From the theory of diffusion we can write: 

D = uj char uj , (9) 

where u c har is the characteristic length in the frequency space. 
For characteristic time in this space we can write: 

tchar = U c har I & = ^ 2 p2 

We can give a physical interpretation for t char for wind-fed pulsars as a 
characteristic time of the momentum transfer: 

For v sw = 1178 km/s we obtain ^ = 245 sec and for v sw = 1281km/s 
we obtain ^ = 200 sec. These velocities are close to v sw obtained using 
estimates of the 'diffusion coefficient'. 

These characteristic time scales are also shown in table 3. 



2.2 Determination of the stellar wind velocity. 

Fluctuating moment $ (see eq.(l)) can be estimated as: 

(Mv t R t \ 



So for D we can write the equation which differs from eq.(6) (Lipunov & 
Popov 1995): 
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Vela X-l 
GX 301-2 
Her X-l 
Cen X-3 



At, sec 

1.5 • 10 4 
1.1 • 10 3 
4 • 10 3 
2 • 10 4 



7 

6.3 • 10 6 
2.9 • 10 6 
6.0 • 10 7 
9.5 • 10 5 



Table 3: 

D, sec -3 

8.7 • 10~ 24 
2 • 10~ 21 
4 • 10- 19 

4.1 • 10~ 17 



v sw , kui/s, 
(v t = 0.1a s ) 

848 

656 



v sw , km/s, 
(v t = a s ) 
1442 
1120 



t chart SeC 

200 
245 
960 
8500 



2 V i / v sw 

where t> t -turbulent velocity, ^-characteristic scale of the turbulence. This 
scale is of order of radius of gravitational capture, Rq- 

where vg = , n8 v , . 

° 10 s cm/ s 

The turbulent velocity is less or equal to the sound speed, a s (in opposite 
case a great bulk of energy will dissipate in the form of shock waves), that's 
why we can write: 

v t = rj -a s , 7] < 1 (14) 

where a s = ((5 RT)/(3 fi)) 1/2 = 1.18-10 6 T 4 1/2 /i _1/2 cm/s and T 4 = T/(10 000AT). 
From eq. (12) we can get: 

D = 4.38 • 10"^^- V^ 2 (j-jy s " 3 . ( 15 ) 

Vio-ro A/f M r J" 

neie iw — 10 i 6g / s , ^45 — i 4 5 s . cm 2 • 

As we see there is a strong dependence of D on ?;. So we can evaluate v (in 

this case it is the stellar wind velocity, v sw ): 

\begin{equation} 
\begin{ array} 
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& v_{sw}=1700\cdot D_{-24}~{-l/7}\dot{M} _{16}~{2/7}\eta ~{2/7}T_4~{l/7} 
\mu ~{-l/7}\\ 

& I_{45}~{-2/7}\left(\frac M{1.5\, M_\odot }\right) "{2/7} \, km/s.\\ 

\end{ array} 

\end{equation} 

Values of v sw are shown in table 3. 

3 Conclusions 

Obtained values of v sw coincide well with the characteristic value of this 
quantity for supergiants of early spectral classes: v sw ~ (600 — 3000) km/s 
(de Jager 1980). We also use well known equation for a terminal velosity of 
a stellar wind which is well confirmed by observations: 

Voo ~ 3 • v esc , (16) 

where v esc = ( 2 ^ M )^ 2 is the escaping velocity on the surface of the star, 
r = _R*.We use this dependence and equation (de Jager 1980): 

v(r) = VU1 - — ) Q (17) 
r 

where a « (0.35 — 0.5). 

For optical component of GX 301-2 we have: M = 35 M Q , R = A3 R Q 
and e = 0.47 (Watson et al 1982). We took r = r min ?a 2i?», because in our 
calculations we used maximum luminosity, i.e. luminosity in periastr. For 
these values we have: v esc ~ 560 km/s and for we have ~ 1680A;m/s 
(for pLeo, Bl lab, in de Jager (1980) we can find v sw = 1580 km/s). 

From eq. (18) we get: v sw (r = 2R*,a = 0.5) ~ 1180 km/s. It coincides 
well with our maximum value « 1150 km/s (see table 3). For Vela X-l there 
are such estimates (Haberl 1991): = 1700 km/ s,r = 1.7 R* and a = 
0.35. We get: v sw ~ 1250 km/s. It also coincides with our estimates of the 
maximum stellar wind velocity: ~ 1450 km/ s. With different assumptions for 
the turbulent velocity (for example v t = |a s ) we can get excellent coincidence 
of our results with the stellar wind velocities estimated above using other 
techniques. 
Acknowledgements. 
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